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FOREWORD 


This report constitutes a part of the work done on the research 
project entitled "Determination of Atmospheric Pollutants from Infrared 
Radiation Measurements." The work was supported by the NASA-Langley 
Research Center through Grant NSG-1282. The grant was monitored by 
Dr. Henry G. Reichle. 
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SUMMARY 

Various molecular band models for infrared radiation are reviewed 
and continuous correlations for the total absorptance of a wide band are 
presented. Different band absorptance correlations are employed in two 
physically realistic problems (radiative transfer in gases with internal 
heat source and heat transfer in laminar flow of absorbing-emitting gases 
between parallel plates) to study their influence on final radiative 
transfer results. 
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1. INTRODUCTION 


The object of this study is to review the molecular band models 
available in the literature [1-14]* and establish their use for atmo- 
spheric applications. Special attention is directed towards obtaining 
the wide band absorptance relations from the basic formulations of the 
narrow band models. Several continuous correlations for the total absorp- 
tance of a wide band are available in the literature [11-15] . These are 
employed in two physically realistic radiative transfer analyses to study 
their effects on the final results of actual radiative transfer process. 

The first problem considered is the problem of radiative energy transfer 
between a gas volume and the surrounding walls. The physical concept of 
this problem can be applied to the radiative transfer analyses of plane- 
parallel atmosphere with non-uniform heating. The second problem considered 
is that of a laminar flow of an absorbing-emitting gas between parallel 
surfaces. Attention is directed, in particular, to carbon monoxide* 
carbon dioxide, methane, and water vapor. 


*Numbers in the brackets indicate references. 



2. BAND ABSORPTANCE AND CORRELATIONS 

The total absorption of a band of overlapping lines strongly depends 
upon. the line intensity, the line half-width, and the spacing between the 
lines. In a particular band, the absorption coefficient varies very rapidly 
with the frequency and, therefore, it becomes very' difficult and time- 
consuming task to evaluate the total band absorptance by numerical integra- 
tion over the actual band contour. Consequently, several approximate band 
models have been proposed [1-14] which represent absorption from an actual 
band with reasonable accuracy. The expressions for absorptance by various 
band models, and several correlations for the totai band absorptance are 
presented in this section. These results are useful in the radiative trans- 
fer analyses presented in Sec. III. 

The absorption within a narrow spectral interval of a vibration rotation 
band can quite accurately be represented by the so-called "narrow band models." 
For a homogeneous path, the total absorptance of a narrow band is given by 

A = 

where K is the mass absorption coefficient, X = p £ is the mass of the 
CO 3 . 

absorbing gas per unit area, £ is the length measured along the direction 

of the path which makes an angle 8 with the vertical, and p_ is the 

density of the absorbing gas. The limits of integration in Eq.(2.1) are 

over the narrow band pass considered. 

Four commonly used narrowband models are Elsasser, Statistical, Random 

* 

Elsasser, and Quasi-Random. The application of a model to a particular case 
depends upon the nature of the absorbing-emitting molecule. Complete discus- 
sions on narrow band models, and expressions for transmittance and integrated 


J ! 

r Aw 


[ 1-exp ( k^X)] 


deo 


( 2 . 1 ) 
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absorptance are available in the literature [1-6,13,14]. 

The total band absorptance of the so-called "wide band models" is 
given by 



[ 1-exp (-KjX)] 


d(to-to Q ) 


t 


( 2 . 2 ) 


where the limits of integration are over the entire band pass and to is 

o 

the wave number at the center of the wide band. 

In actual radiative transfer analyses, the quantity of frequent interest 
is the derivative of Eqs.(2.1) and (2.2). With reference to Eq.(2.2), this 
can be expressed as 


dA/dX * A' (< W ,X) 



exp (-K^X) ] d(to-to Q ) 


(2.3a) 


For multiband gases, Eq. (2.3a) can be written as 


a' - E A’0c wi ,X) = E / exp (-K U .X) ] d(» 

i=l 1=1 Ato^ 


(2.3b) 


where Ato^ is the spectral range of the ith wide band and n represents 

the number of important bands of the absorbing-emitting gas under consideration. 


2 • 1 Exponential Wide Band Absorptance from Narrow Band Models 

Detailed discussions on the wide band models are available in the 
literature [7-15] . The important formulations pertaining only to the exponen- 
tial wide band models are presented here. The relations for integrated (total) 
band absorptance for these wioe bands are obtained from the absorptance 
formulations of. narrow band models by applying the following relation for the 
variation of line intensity [10,13-15], 

iWnOBUCtBtLITY OF Till 
PAGE IS POOP, 
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S^/d = (S/A o ) exp {[-b o |w-a) o |l/A o } , (2.4) 

where is the intensity of the jth spectral line, d is the line 

spacing, S is the integrated intensity of a wide band, A q is the band 

width parameter, and b -2 for a symmetrical band and b =1 for bands 

o o 

with upper and lower wave number heads at u^. The total absorptance of a 
exponential wide band is obtained from the absorptance of narrow band models 
by employing the relation 

A(u,3) - A(u,3)/A n = f [L, (u,3)] d(o»-to) ) , (2.5) 

•'wide ° 

band 

where u = SX/A is the nondimens ional path length, 3 = 27iy /d is the 
line structure parameter, y L is the Lorentz line hald-width, and A^(u,3) 
represents the mean absorptance of a narrow band. 

By employing the Elsasser narrow band absorptance relation and Eq.(2.4), 
the expression for the exponential wide band absorptance is obtained as (13,14] 


A(u,3) = Y + (1/TT) 



[An iJj + E^(ij/)] dz 


( 2 . 6 ) 


where ip - u sinh 3/(cosh 3 - cos z) , y = 0.5772156 is the Euler’s constant, 
and E^OJj) is the exponential integral of the first order. Analytic solu- 
tion of Eq.(2.6) can be obtained in a series form as [13] 


A(u , 3) = £ {-(A) n [SUM(mn)]/[n(B+l) n n! (n-1)!]} , 
n=l 


where „ 

00 

SUM(mn) = [ (n-fm-1) ! (2m-l) ! C m ] / [2 m (m! ) 2 ] , 

m=o 


(2.7) 


A ** -u tanh 3» B = 1/cosh 3 » C = 2/(1 + cosh 3) = 2B/(B+1) . J 
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It can be shown that the series in Eq.(2.7) converges rapidly. When the 
weak line approximation for the Elsasser model is valid (i.e., $ is large), 
then Eq.(2.6) reduces to [13,14] 

A(u) = y + in (u) + E^(u) . (2.8) 

In the linear and logarithmic limits, Eqs.(2.6) and (2.7) reduce to 


A ■ u , (u « 1 , 3 > 1) 

A * y + £n(u) , (u » 1 , 00 > 1) • 


(2.9) 


It can be shown that Eq.(2.6) reduces to the correct limiting form in the 
square-root limit [14J. 

By employing the uniform statistical narrow band absorptance relation 
and Eq.(2.4), the expression for the exponential wide band absorptance is 
obtained as [13,14] 


where 


A(u 


,6) - f 1 U- 


exp [-3 L(n)]} £ 1 d£ , 1 


n = u^/8 , ? = expt-bJw-wJ/A ] . 


( 2 . 10 ) 


In the limit of large 3 , this reduces to Eq.(2.8), and the linear and 

logarithmic limits are given by Eq.(2.9). This also reduces to the correct 

limiting form in the square-root limit [14]. 

By employing the general statistical narrow band absorptance relation 

► 

and Eq. (2.4), a third expression for the exponential wide band absorptance 
is obtained as [13-15] 


pjSRODuemiLmr op 

onr-iMl Mas is pot» 
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A(u,8) - f\[l-exp<- t )]/[5 2 + <2p u ) 2 l 1/2 > d5 
•'o 

+ f { [1-exp (-pt)]/£}d£ , 
where p = [l+(u^/t>] - l/[l+(u£/t)) , 

p u = {(t/u)[l+(t/u)]r 1/2 , 


( 2 . 11 ) 


S = p/p u , t = 0/2 * ir y L /d . 


In the linear and logarithmic limits, this reduces to the expressions given 
by Eq.(2.9). 

A fourth expression for the exponential wide band absorptance is 
obtained from the narrow random Elsasser model as [13,14] 


A(u,$) = £n{[G 2 (u,$) + u sinh 3 + cosh 3]/(sinh 3 + cosh 3 )) » (2.12) 

where 

C^u^) = [(u 2 + 1) sinh 2 3 + 2u sinh 3 cosh 33^ 2 . 

* 

It may easily be shown that Eq.(2.12) satisfies the linear, square-root 
and logarithmic limits for appropriate u values [14]. 

The solutions of the four wide band absorptance listed above are 
termed as the exact solutions. These are illustrated in Fig. 2,1 for three 
different values of the line structure parameter 3 • The agreement between 
the results of various models are found to be within 25%. Since the absorp- 
tion at small path lengths is a function solely of the total intensity of 
absorbing lines, the results by all models agree well in this region. At 







larger path lengths, however, the Elsasser model predicts higher absorption 
than the general or uniform statistical model. This is because there is 
always more overlapping of the spectral lines in the statistical models 
than in the regular Elsasser band. The main reason for presenting these 
results here is because various correlations based on these results are 
employed in actual radiative transfer analyses in Sec. III. 


2.2 Band Absorptance Correlations 

Several continuous correlations for the total absorptance of a wide band, 
which are valid over different values of path length and line structure para- 
meter, are available in the literature. These are discussed, in detail, in 
[11-15] and are presented here in the sequence they become available in the 
literature. Most of these correlations are developed to satisfy at least some 
of the limiting conditions (nonoverlapping line, linear, weak line, and strong 
line approximations, and square-root, large pressure, and large path length 
limits) for the total band absorptance [12,14]. Some of the correlations even 
have experimental justifications [8,11]. The results of these correlations 
are compared in Fig. 2.2 for t = 0.01 and 1. Figure 2.3 shows the errors 
encountered in using the various correlations as compared with the exact solu- 
tion obtained by the general statistical model, Eq.(2.11). For comparison of 
other results, one should refer to [13,14]. 


(a) Tien and Lowder [11]: 

A(u,0) = fn(uf(t) { (u+2) / [u+2f ( t) ) ] } +1) , 

where f(t) = 2.94 [1-exp (-2.60 t)jj , t = 0/2 . 


(2.13) 


This correlation does not reduce to the correct limiting form in the 
square-root limit [12] , and its use should be made for 0 > 0.1 . 


(b) Goody and Belton [16] : 

A(u,0) = 2 in {1 + u/[4 + (im/4t)] 1/2 ) 



m 10 

\ o, ■ PAGE IS PGUP: 
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The use o£ this correlation, should be made for relatively small 3 values 

[ 12 ]. 

(c) Tien and Ling [17]: 

A(u) = sinh l (u) . (2.15) 

This relation is valid only for the limit of large 3 • 


(d) Cess and Tiwari [12-14] : 

■ . ' . : . : i 

A(u,3) = 2 £n(l+u/{2+[u(l+l/3)] 1/2 }) , 

where _ 

3 = 4t/iT = 23 /tt . 

•: • I 

The use of this correlation is justified at relatively high pressures to 
gases whose spectral behavior can be described by the general statistical 
model. By slightly modifying Eq.(2.16) another form of the wide band 
absorptance is obtained as [13,14] 


(2.16) 


A(u,3) = 2 2n(l+u/{2+[u(c+ir/23)] 1/2 }) 


where 


0.1 » 

8 

< 

1 

and 

all 

u values 


c =* " 

0.1 , 

e 

> 

1 

and 

u £ 

1 



0.25 , 

8 

> 

1 

and 

u > 

1 


(2.17) 


Equations (2.16) and (2.17) reduce to all the limiting forms [12]. The 
results of these correlations are compared with other correlations in 
Figs. 2.2 and 2.3 . 

(e) Edwards and Balakrishnan [10]: 


A(u) = £n(u) + E^u) + y + — - E 3 (u) . 


(2.18) 


This correlation is valid for large 3 and for large values of path lengths 
(see Figs. 2.2 and 2.3). 
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(f) Felske and Tien [15]: 


A(u,8) = 2 E x (tp u ) + E x (P u /2 - E 1 [(p u /2)(l+2t)] 


+ £n[ (tp u ) 2 / (l+2t) ] + 2y 


(2.19) 


This correlation is valid for the entire range of the governing parameters 
(Fig. 2.2). The comparison of results of this correlation with the exact 
solution of Eq.(2.11) indicates excellent agreement (Fig. 2.3). This, however, 
would be expected because the correlation was obtained from Eq. (2.11). 

(g) Tiwari and Batki .[13,14]. The form of the absorptance given by Eq. (2.8) 
can be treated as another correlation for the total band absorptance. The use 
of this correlation is justified at all path lengths for t - (8/2) 1 (Fig. 2.2) 

A comparison of results presented in Figs. 2.1 and 2.2 indicates that Tien 
and Lowder's correlation follows the general trend of the Elsasser model while 
Cess and Tiwari' s correlation follows the general statistical model. This fact 
was noted clearly from several comparative results presented in [14]. It was 
indicated in [14] that results of various correlations could be in error up to 
40% when compared with the exact solutions based on different band models. Tien 
and Lowder's correlation was found to give the least error when compared with 
the exact solution based on the Elsasser model for t = 0.1, 1, and 10. 

From the comparsion of results presented in Figs. 2.2 and 2.3 (and from 
results of references [13,14]), it may be concluded that the correlation of 
Felske and Tien provides fairly accurate results for all pressures and path 
lengths. At relatively high pressures, however, the correlation given by 
Eq.(2.8) provides a uniformly better approximation for the total band absorp- 
tance. While Felske and Tien's correlation is suitable for radiative transfer 
analyses in gases exhibiting statistical spectral behavior, a better correla- 
tion is needed for gases whose spectral characteristics warrant use of the 
Elsasser model. One such relation is given by Eq.(2.7) which is the series 

form solution of the Elsasser total band absorptance. 


— LARGE PATH LENGTH LIMIT 

— : LINEAR LIMIT 

t a 0.01 


t 3 1.0 


t«/5/2=7ry/d 
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Fig. 2.2 Comparison of results of band absorptance correlations for t * 0.01 and 0.1 





u» Spy /A 0 


Fig. 2.3 Errors in the band absorptance correlations when compared with the 
exact solution (General Statistical Model) for t = 0.01 and 0.1. 


3 . RADIATIVE TRANSFER ANALYSES 

The band absorptance correlations discussed in the previous section 
are employed in this section to two illustrative radiative transfer problems. 
The physical model and coordinate systems for both problems are shown in 
Fig. 3.1 . 

The expression for the total radiative flux may, in general, be given 
by [12,18] 


*R 


w 

‘ / ’ro. d “ - e i ' e ; 


+ 1 / t e u (z > - < T i> ] / K „i ®'pt-iv y ' z)1 d “ iz 


o o 
L 

y © 


Aw 


- I / l e u <z > - e w (T ?. >] / k u ex P [ -| K o, (z_!r)I dz - 

v o o Au) 


(3 


where e = aT 4 ', with a denoting the Stefan-Boltzmann constant, and Au> 
indicates integration over the single band. The derivation of Eq.(3.1) 
assumes the existence of a local thermodynamic equilibrium and employes the 
exponential kernal approximation. The Planck function, e^(T) , is considered 
a slowly varying function of wave number over the narrow band and its value 
is evaluated at the band center. The primary motivation for employing the 
exponential kernal approximation is that it allows the kernal function in 
Eq.(3.1) to be expressed in terms of the derivative of the total band 
absorptance which is given by Eq.(2.3). 

A combination of Eqs.(2.3) and (3.1) results in 

q R CO =* e x -e 2 + f A^C / [e^ (?') - e^ df-' 


- { [ \ (V) ’ V CT 2 )IX,[ ! dD , 


(3 


Eligibility of tils 
L PAGE IS POOR 



where £ « y/L , u = S(T ) p L/A (T_), A(u,£) » A(u,£)/A . 

O A O X O 

In this equation. A* (u) denotes the derivative of A(u) with respect 
to u and u q represents the nondimens ional pressure path length. 

Equation (3.2) possesses two limiting forms (optically thin and large 
path length limits) which are discussed in detail in [12,19,20], It 
should be pointed out here that while the large path length limit may 
depend upon a particular band model employed in the radiative flux equation, 
the optically thin limit is completely independent of the band model. 

In order to study the effect of different band absorptance correla- 

...... i 

tions on the final results of actual radiative transfer processes, two 
illustrative cases (extensively studied in the literature) are treated here. 
The first is the problem of radiative energy transfer between a gas volume 
(where there is a uniform heat source) and the surrounding walls. The phys- 
ical concept of this problem can be applied to the radiative transfer analyses 
of plane-parallel atmosphere with non-uniform heating. The second problem 
considered is the problem of liminar flow of absorbing-emitting gases between 
parallel surfaces. This represents a realistic situation of convective- 
radiative energy transfer in the atmosphere and in other engineering problems. 
For mathematical simplicity, the bounding surfaces for both problems are 
considered to be black. 

3.1 Radiative Transfer in Gases with Internal Heat Source 

The physical model and the coordinate system for this case is shown in 
Fig. 3.1a , The two parallel black walls are considered to be at the same 
constant temperature T^ . The gas is assumed to have a uniform heat source 
(ot sink) of strength Q per unit volume. The radiative transfer is the 
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sole mechanism of energy transfer within the gas. The local temperature 
distribution is thus a consequence of the uniform heat source, adding 
energy to the gas which in turn is transferred through the gas to the bounding 
surfaces by radiative transfer. 

For this special case the conservation of energy yields 

dq R /dy = Q (3.3) 

Since the problem is symmetric, an integration of Eq.(3.3) gives 
q «= (1/2) Q L (2£ - 1) and after combining this with Eq.(3.2) results in 

« - f = f{ j [|u o «-5’)] ds* 

o 

- / <kd A &U(V<)1 dt'} , 

£* ^ o 

where <KS) = [« w (T) - e w (^>1 / (Q/PS (T )) . 

o o 

In this equation <p(Z) represents the temperature profile within the gas 
in terms of the Planck function. The equation is written for a single band 
gas but it can easily be extended to the case involving multi-band gases [12]. 

In the optically thin limit (i.e., for u q « 1) , A(u) = u , and 
A (u) =» 1 such that Eq.(3.4) yields the result [12,20] 

- 1/3 . (3.5) 

For most band models, the band absorptance in the large path length 

limit (i.e., for u » 1) is given by the relation A(u) = Jin (u) such 

o 

that A (u) = 1/u . In this limit, therefore, Eq.(3.4) reduces to [12,19] 

z - J 1 <P(V) d V/<£-V) , 

o 



(3.6) 



for which the solution is found to be 


KO - a/ir)E«l - 01 1/2 . (3.7) 

As discussed in [12,19J, Eq.(3.7) yields the result that the gas temperature 
at a surface is equal to the surface temperature. This absence of a temp- 
erature slip is characteristic of optically thick radiation. 

The first information needed for numerical solution of Eq.(3.4) is 
the derivative of the band absorptance correlations employed in the analysis. 
Numerical solutions of (3.4) are obtained by using the method of undetermined 
parameters. In this method a polynomial solution for <J> is assumed, and the 
constants are evaluated by satisfying the integral equations at equally 
spaced locations. If a quadratic solution for <f> is assumed then 
<KO * a Q + + a 2 » ant * by using the symmetry condition <f> ' (£= 1/2) = 0 , 

one obtains 

<KS) - a Q + a L (S - S 2 > . (3.8) 

The constants a^ and a^ are evaluated by satisfying Eq.(3.4) at £ ** 0 
and 1/4 . If quartic solution for <J> is assumed then by using the symmetry 
conditions one obtains 

<KO - a x - a 2 a - £ 2 ) - a 3 (| £ - £ 3 +j £ 4 ) . (3.9) 

The constants a^, a 2 , and a^ are obtained by satisfying Eq.(3.4) at 
£ = O', 1/6, and 1/3 . The choice of locations for satisfying the integral 
equation to obtain the constants is consistent with the symmetry conditions. 

By employing the various band absorptance correlations discussed in the 
previous section, numerical solutions of Eq.(3.4) were obtained for different 
values of the line structure parameter t =8/2 . The results are illus- 
trated in Figs.’ 3.2 - 3.4 along with the limiting solutions. It should be 
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pointed out here that these results apply to any situation for which radiative 
transfer within the gas is the result of a single band. 

It is obvious from Figs. 3. 2-3. 4 that the band absorptance results 
of all the correlations approach the optically thin limit for small u Q , 
the influence of the line structure parameter is maximum for intermediate 
values of u Q , and in the large path length limit the solutions become 
independent of g . The one exception to this is that the results of Goody 
and Belton's correlation do not approach the correct logarithmic limit for 
large g . The reasons for this are discussed in [12,14] where it was pointed 
out that the use of the Goody and Belton's correlation should be restricted 
to relatively small values of g . 

For a particular value of g , the results of different correlations 

approach the linear and logarithmic limits for different u Q values. For 

g = 0.02, for example, the results of all correlations (except Tien and 

Lowder's) almost are identical for u > 1. This corresponds to the range of 

o — 

the square-root limit where three separate conditions (g « 1, u/g » 1, and 
g « 1) must be satisfied [12,14]. The square-root limit is not satisfied 
by the Tien and Lowder's correlation. Although the linear limit is independent 
of any spectral model, it is approached by different correlations for different 
u Q values. 

Maximum differences between the results of various correlations occur for 
the intermediate values of u Q (0.1 < u c < 10) and for 0.1 < g < 1 . For large 
g values (i.e., for high pressures), the line structure of a gas is smeared 
out and differences between the results of various correlations become small 
(see Fig. 3.4). This situation corresponds to the weak line approximation 
[14], and some consideration of this usually is given in developing a particular 
correlation. 
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In general (and for the intermediate path lengths in particular) , the 
temperature difference between the centerline of the physical system and 
the wall is found to be lowest for the Tien and Lowder's correlation and 
highest for the Cess and Tiwari's correlations. For the most part, the 
results of other correlations fall between these two extreme values. The 
physical reasoning behind this is consistent with the results of the total 
band absorptance illustrated in Fig. 2.2. If the absorption by a gas is high, 
then, to maintain the local thermodynamic equilibrium, the emission by the 
gas will be higher. This will result in high radiative energy transfer from 
the gas and consequently will lower the center line temperature. The trend 
exhibited by the results of Fig. 2.2 and Figs. (3.2 and 3.3) is, therefore, 
mutually consistent. Thus, the difference in the centerline temperatures, as 
obtained by using the different correlations, are of the same order as the 
difference in the band absorptance results of various correlations. 

As mentioned in Sec. 2 (see Fig. 2.1), the Elsasser theory always pre- 
dicts higher absorption than the general statistical model. From this and the 
results of Figs. 3.2 - 3.4, it may be concluded that for 3 > 0.02, use of the 
Tien and Lowder's correlation is justified for gases whose spectral behavior 
can be described by the regular Elsasser model. On the other hand, use of 
the Cuss and Tiwari's correlation could be made to gases with bands of highly 
overlapping lines. Since, in a vibration-rotation band, the lines are dis- 
tributed neither regularly nor randomly, care should be taken in applying a 
particular correlation for a specific case. 
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Fig. 3.2 Comparison of results for a single band gas as obtained by using the various band 
, absorptance correlations. 
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Fig. 3.3 Comparison of results for a single band gas obtained by using the various band 
absorptance correlations for t = 1. 




HCIBILITY OF THE 
;L PAGE IS POOR 



Fig. 3. A Comparison of results for a single band gas obtained by using the various band 
absorptance correlations for t = 10. 


K> 

u> 




24 


3.2 Heat Transfer in Laminar Flow of Abosrbing-Emitting Gases 

Between Parallel Plates 

The second problem considered is the problem of heat transfer in 
laminar, incompressible, constant properties, fully developed flow of 
absorbing-emitting gases between parallel plates. The physical model and 
coordinate system for this case is shown in Fig. 3.1b . The conditions of 
uniform surface heat flux for each plate is assumed such that the tempera- 
ture of the plates, T^ , varies in the axial direction. Fully developed 
heat transfer is considered, and axial conduction and radiation is assumed 
to be negligible as compared with the normal components. Consistent with 
the constant properties flow, the absorption coefficient is taken to be 
independent of temperature and radiation can be linearized. Extensive treat- 
ment; of this problem is. available in the literature [8,12,21-27] . The sole 
motivation for studying the problem here is to investigate the influence of 
different band absorptance correlations on the radiative transfer capability 
of a gas in a more realistic situation. 

Within the confines of foregoing assumptions, the energy equation for 
the present problem can be written as 

3T d Zr £ i 

V X 35T ’ Q W ' pS~ ( V 3y) • < 3 - 10 > 

p 

where the parabolic velocity profile is described by the equation 

V x = 6 v m Ky/L) - (y/L) 2 ] . < 3 - 11 > 

The mean velocity v is given by the expression 

V m = I" / v x dA c • (3.12) 

G A 

C 

where A is the cross-sectional area normal to the direction of flow, 
c 

For a uniform wall heat flux and fully developed heat transfer £T/9x 
is a constant and is given by the expression 

- 2 o q w /(v m LX). (3.13) 
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where X is the thermal conductivity of the gas. 


A combination of Eqs.(3.10) through (3.13) results in 


12(E-£ 2 ) - ^ 

4 ; dF d? » 


V 


where 


(3.14) 


C = y/L , 0 - (T-Tj)/ (q^L/X) . 

Upon integrating this equation once, and by noting that d0/d£ = 0 and 


q * 0 at £ = 1/2 , one finds 
R 


de 

d£ 


- 2 (3? 2 - 2C 3 ) + 1 = q R /q w 


(3.15) 


By combining Eqs.(3.2) and (3.15) the governing integro-dif f erential 
equation for linearized radiation in multi-band gases can be written as [12] 


d0 

dS 


- 2(3£ 2 - 2£ 3 ) + 1 =* 


lx t H i V { f e( ?'> \ «• 

X— 1 o 

- / 0(C) sj [|u 0i «'-O] dS’> , 


n 


where 


h = Ai /A oi , U Qi = S ± (T 1 ) PL/A^dj) 


H i ‘ A <,i< d< W dT > Il • H - 2 «i • 


In this equation S^(T^) is intensity and A q ^(T^) is t ^ ie k an< * 
width parameter for the ith band. The equation describes the temperature 
profile within the gas for which the boundary condition can be written as 
0 ( 0 ) - 0 . 

For flow problems, the quantity of primary interest is the bulk 
temperature of the gas, which is defined as 


(3.16) 
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mA 


/ 


T v dA 
X c 


JL / 

v L J 
m o 


Tv dy 
x 7 


(3.17) 


By employing Eqs.(3.11) and (3.12) and the nondimensional quantities 
defined earlier, the bulk temperature can be expressed in a dimensionless 
form as 

-1 

0 b = (T b - T.^) / (q^ L/X) = 6 J 0(5) (5 - £ 2 ) d£ . (3.18) 

The heat transfer q^, is given by the expression, q w = h c (T^ - T b ) , 

where h is the convective heat transfer coefficient (btu/hr-sq ft*-°F). 
c 

In general, the heat transfer results are expressed in terms of the Nusselt 
number Nu, a dimensionless quantity defined as, Nu = h c D b /X . Here, D b 
represents the hydraulic diameter, and for parallel plate geometry it equals 
twice the plates separation, i.e., D b = 2L „ Upon eliminating the convective 
heat transfer coefficient, h c , from the expressions for q^ and Nu , a 
relation between the Nusselt number and the bulk temperature is obtained as 

Nu - 21 - T b ) = -2/0 b . (3.19) 

The general solution of Eq.(3.16) is obtained by numerical procedures. 

Before discussing the methods of obtaining these solutions, it Would be 
convenient to discuss the limiting forms of the general equation. 

For the case of negligible radiation, the divergence of radiative flux 
is taken to be zero and Eq.(3.14) yields 

- 0 ( 5 ) = 5 - 25 3 + 5 4 • .(3.20) 

By substituting this into Eq.(3.18), the bulk temperature for this case of 
no radiation is obtained as 

(3.21) 


0 = - 17/70 

o 



In the optically thin limit Eq.(3.l6) reduces to 


d6 


- 2(3? 2 - 25 3 + 1) 


n 


+ix t h v { r d «'- j -1 ««’) dc-j . 

i=l o 5 


(3.22) 


After differentiating once, Eq.(3.22) can be expressed in the following form 


d 2 0 

dp 


• 3 N 9(5) = 12(5~5 2 ) , 


(3.23) 


for which the boundary conditions are 9(0) - 0 , and 0'(l/2) = 0 . 

The quantity N in this equation is referred to as the optically thin 
radiation-conduction interaction parameter and is given by the relation 


n 


N = <p L 2 /X) £ S (T ) (d e(0l /dT) T 
i=l 1 


(3.24) 


Equation (3.23) possesses an elementary solution, and the result expressed 
in terms of the bulk temperature is found to be 

0 b = [1/ (3N) 3 ] { (576/v' / 3N) (NEXP) - 21.6 N 2 + 72 N - 288} , 

NEXP = [1 - exp(-/3N)]/[l + exp(-/3N)] . 


where 


(3.25) 


The results of this equation (0^ versus N) are illustrated in Fig. 3.5 
along with other limiting solutions. 

In the large path length limit the band absorptance for each band of 
importance is given by the relation A^(u^) = n(u^) such that A ^(u^) = 
1/u^ . Upon substituting this in Eq.(3.16) there is obtained 


where 


d0 

d5 


-x 

- 2(3? 2 - 25 3 ) = M J 0(5') dZ'/tt-V) 


M “ HL/X = (L/X) £ -(L/X) £ A oi (de a,i /dT) T 
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(3.26) 


28 


The dimensionless parameter M constitutes the radiation-conduction inter- 
action parameter for the large path length limit. Equation (3.26) does 
not appear to possess a closed form solution. A numerical solution has 
thus been obtained, and the results for bulk temperature (0^ versus M) 
are illustrated in Fig. 3.5 . 

Numerical solutions of Eqs.(3.16) and (3.26) are obtained by the 
method of undetermined parameters. For this case, a polynomial solution 
for 9(5) is assumed as 

6(5) - a Q + a^5 + + a $ ^ + ^ ’ (3*27) 

After employing the conditions 9(0) — 0 , 0'(l/2) = 0 , and 9'(1) ™ 

- 0'(O) , Eq. (3.27) becomes 

0(5) = a 1 (5 - 25 3 + 5 4 ) ■+ a 2 (5 2 - 25 3 + 5") • (3-28) 

The constants a 1 and a 2 are obtained by satisfying the governing 
integro-differential equation at two locations 5 “ 0 and 5 “ 1/A . 

A combination of Eqs.(3.18) and (3.28) results in 

e b " 77 (17a l + 3a 2 5 • (3 - 29) 

Thus, with a, and a known, the bulk temperature (or Nusselt number) 

1 2 

is obtained from Eq. (3.29) . 


)/ (qJ-/A.) 


..4.. 


(N = PL^/A) or (M=HL/A) 


Fig. 3.5 Limiting solutions of the flow problem. The abscissa for optically thin limit 
is N and for large path length limit is M. 
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A complete discussion (and physical interpretations) of the various 
parameters entering into the present problem is given in [12,22,23). Numeri- 
cal solutions of Eq. (3.16) were obtained in terms of the bulk temperature 
(see Eq. (3.29)), and these are presented here in terms of the dimensional 
quantities L and P. Specific results were obtained for CO (fundamental + 

1st overtone bands), C0 2 (15y, 4.3y, and 2.7y bands), H 2 0 (rotational, 6.3y, 
2.7y, 1.87y and 1.38y bands), and CH 4 (7.6y and 3.3y bands) for which spectral 
information were obtained from [8,11,23). 

For CO, C0 2 , H 2 0, and CH 4 , the results obtained by employing the Tien 
and Lowder’s correlation for band absorptance are illustrated in Figs. 3.6 
through 3.10. The limiting value of 0^ = -0.243 corresponds to negligible 
radiation, and the effect of radiation increases with increasing plate spacing. 
As would be expected, radiative transfer is more pronounced for higher pressures 
and wall temperatures. Also shown in Figs. 3.6 - 3.9 are the limiting solutions 
for large u Q (LLU) . It is seen that, for a given wall temperature, the large 
u Q limit can be obtained either by going to large values of L or to high pres- 
sures. These results also indicate that, for a particular wall temperature, 
the large u Q limit for C0 2 is achieved at a relatively lower pressure than for 
other gases; As a matter of fact, for most practical purposes involving C0 2 
at room temperature, the results for one atmosphere can be regarded as the 
results for the large u Q limit. A comparison of the results for the four 
gases is shown in Fig. 3.10 for a pressure of one atmosphere and a wall temp- 
erature of 1,000 °K. As discussed in [12,19,20], the relative order of the 

• :■ : I ' ■' ' • • : ; 

four curves, for small values of L is characteristic of the interaction 
parameter for optically thin radiation, and for large L is is characteristic 
of the interaction parameter for the large u Q limit. 
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Fig. 3.7 Results for CO 2 obtained by using the Tien and Lowder's correlation. 
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Bulk temperature results for CO (fundamental band) and C0 2 (15y, 4.3y, 
and 2.7y bands) as obtained by employing the various correlations for band 
absorptance, are illustrated in Figs. 3.11 through 3.14. Results for CO are 
illustrated in Figs. 3.11 and 3.12 for wall temperatures of 500 °K and 1,000 °K 
respectively. It is evident from these figures that, except for the results 
of Tien and Lowder's correlation, results of other correlations differ from 
each other by less than 6% for all pressures and path lengths. For P * 0.1 
and 1 atm, results differ by not more than 3%. The largest difference of about 
6% occurs for P = 10.0 atm and T x = 1,000 °K. From a close observation of all 
the results presented in Figs. 3.11 and 3.12, it may be concluded that, for 
low to moderate pressures (say up to 5 atm), any one of the correlations 
(No. 2,3,5, or 6) could be employed in radiative transfer analyses. At high 
pressures, however, use of correlations 5 or 7 is recommended. 

From Ref. [22] and the results presented in Figs. 3.6, it is noted that 
for CO the limit of large u Q is approached at P - 10 atm for Tj = 500 °K and 
at 100 atm for T : = 1,000 °K. This trend is also evident, in general, from 
the results of Figs. 3.11 and 3.12. The results of Tien and Lowder's correla- 
tion, however, follow this trend more closely than any other result. As such, 
use of the Tien and Lowder's correlation is justified for radiative transfer 
analyses involving gases like CO (i.e., diatomic gases with single fundamental 
band having uniform distribution of spectral lines) at moderate and high pres- 
sures . 

For C0 2 , results of different correlations are illustrated in Figs. 3.13 
and 3.14 for P = 0.01, 0.1, 1, and 10 atm, and for T J = 500 °K and 1,000 °K 
respectively. From Ref. [23] and the results presented in Fig. 3.7, it may 
be noted that for C0 2 the limit of large u Q (LLU) is approached at 2 atm for 



T x = 300 °K, at about 4 atm for T : = 500 °K, and at about 10 atm for T 1 =1,000 °K. 
Thus, results for 10 atm in Figs. 3.13 and 3.14 essentially are LLU results. 

For clarity, results of P = 1 and 10 atm are not plotted on the same graph. 

As was the case with CO, the results of all correlations (except Tien 
and Lowder’s) almost are identical for CO 2 also for P = 0.01 and 0.1 atm. 

This, however, would be expected because the low pressure (small 3) situation 
corresponds to the case of square-root limit and most correlations are devel- 
oped to satisfy this limit. It was pointed out earlier and in [12,14] that 
the square-root limit is not satisfied by the Tien and Lowder's correlation. 

At low pressures, therefore, use of the Tien and Lowder's correlation cer- 
tainly is not justified. Other results of C0 2 (shown in Figs. 3.13 and 3.14) 
follow the same general trend as for CO in Figs. 3.11 and 3.12. The maximum 
difference between the results of different correlations is about 6% for 
P = 1 atm and Tj = 1,000 °K. For the most part, results of correlations 3, 

5, 6, and 7 are identical for P - 10 atm. This again would be expected 
because for C0 2 , the LLU is approached at relatively lower pressures and most 
correlations are developed to satisfy the logarithmic limit. For gases like 
CO 2 , therefore, use of any one of the correlations 2, 3, 5, and 6 is recommended 
at low and moderate pressures, and of 3, 5, 6, and 7 at high pressures. Use of 
the correlations 2, 3, 6, and 7, in a particular radiative transfer analysis, 

t 

provides a greater mathematical flexibility and simplicity. 
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Fig. 3.11 Results for CO (fundamental band) with Ti = 500 K. 
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Fig. 3.13b Results for CO 2 (three bands) with Tj = 500 °K. 
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Fig. 3.14b Results for CO 2 (three bands) with Ti * 1,0G0 °K 





CONCLUSIONS 


In this study, use of several continuous correlations for total band 
absorptance were made to two problems to investigate their influence on the 
final results of actual radiative processes. -For the case of radiative 
transfer in a gas with internal heat source, it was found that actual 
centerline temperature results obtained by using the different correlations 
follow the same general trend as the results of total band absorptance by 
these correlations. From these results, it may be concluded that use of the 
Tien and Lowder's correlation should be avoided at lower pressures, but its 
use is justified (at moderate and high pressures) to gases whose spectral 
behavior can be described by the regular Elsasser band model. For all 
pressures and path length conditions, use of the Cess and Tiwari's corre- 
lations could be made to gases with bands of highly overlapping lines. In a 
more realistic problem involving flow of an absorbing-emitting gas, results 
of different correlations (except the Tien and Lowder’s correlation) differ 
from each other by less than 6% for all pressures and path lengths. Use of 
Tien and Lowder's correlations is justified for gases like CO at moderate 
and high pressures. For gases like CO 2 , use of any other correlation is 
recommended. While Felske and Tien's correlation is useful for all pressures 
and path lengths to gases having random band structure, Tiwari and Batki's 
simple correlation could be employed to gases with regular or random band 
structure but for P > 1.0 atm. 
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